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The BaBar measurements of the ratios 7?.(D'*h = deviate from the standard model 

^ ' B(S->i3l*Vi'(») 

expectation, while new results on the purely leptonic B — ^ rPr mode show a better consistency 

with the standard model, within the uncertainties. In a now physics scenario, one possibility to 
accomodate these two experimental facts consists in considering an additional tensor operator in 
the effective weak hamiltonian. Wc study the effects of such an operator in a set of obscrvables, in 
semileptonic B — ^ Z)'*^ modes as well as in semileptonic B and B^ decays to excited positive parity 
charmed mesons. 
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I. INTRODUCTION 

The BaBar measurements of the rates of and semileptonic decays into I?^*' and a r lepton seem to indicate a 
significant deviation from the standard model (SM) expectation. The experimental results concern the B -J> 
decay widths normalized to the widths of the corresponding modes having a light i = e, fj, lepton in the final state [T]: 

B(B" D^T^ 9 ) B(B^ ^ D*^T^ V ") 

n\D) . ;(^o^^.,-,:; - 0-469 ± 0.084 ± 0.053 , 7^"p*) ^ - 0.355 ± 0.039 ± 0.021 

(1) 

(the first and second error are the statistic and systematic uncertainty, respectively). The measurements have been 
estimated to deviate at the global level of 3.4(t with respect to SM predictions [HH]. Therefore, there is the possibility 
that semileptonic processes involving heavy quarks and the r lepton are unveiling the efi^ects of particles with large 
couplings to the heavier fermions, as it is natural for charged scalars which could contribute to the tree-level h — d.v 
transitions [2l-[9]. 

Before the observation of these possible hints of new physics (NP) in semileptonic 6 — > c decays, the first experimental 
analyses of the purely leptonic B~ t~Vt mode also reported an excess of events. In SM the B{B~ t~Vt^ 
branching fraction is given by 



B{B- ^ r-D,) = "^J - [I a iKbl^r^^ , (2) 

neglecting a tiny electromagnetic radiative correction. Using the lattice QCD average for the B decay constant 
fs ~ (190.6 ± 4.7) MeV quoted in [10], and varying the Cabibbo-Kobayashi-Maskawa (CKM) matrix element \Vub\ 
in the range determined from inclusive and exclusive B decays: iT^bl = 0.0035 ± 0.0005, the prediction follows: 
B{B- T'Vr) ^ (0.79 ± 0.23) x 10""^, in agreement with the outcome of CKM matrix fits [TTl [H]. This value 
is smaller by about a factor of 2 than the experimental results reported in [TB1 - [TC| and compiled in [XT': B{B~ — 
T~Vr) = (1.68 ± 0.31) X lO^"'. However, new Belle [18] and BaBar [19] measurements, obtained using the hadronic 
tagging method, 

B{B- -> T-9r) = (0.72i°|^ ± 0.11) X 10""* (Belle) 

B{B- T-Pr) = {iMtoil ± 0.24) X lO"'' (BaBar) (3) 

are more consistent with SM, and draw the average B{B~ t^Vt) to a smaller value: B{B^ T^Vr) ~ (1.12 ± 
0.22) X lO"'*, after the combination with the semileptonic tagging method results, see fig[l] 
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FIG. 1. Experimental results for B{B~ — >■ T~Vr) [141 1161 [THl 119) together with the SM expectation corresponding to \Vub\ = 
0.0035 ± 0.0005 (vertical band). 

The different trend of the measurements involving r in i? leptonic and semileptonic decay modes poses two questions. 
The first one concerns the level of accuracy of the SM predictions for the ratios in ([I]). The second one is which kind 
of new physics effects, if any, could modify the ratios ([T]) without affecting the purely leptonic mode. Indeed, several 
analyses devoted to try to explain the anomalies in _B — > D^*^TPr within new physics scenarios have considered as 
possible candidates models with new scalars having couplings to leptons proportional to the lepton mass, to guarantee 
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the enhancement of the r modes. This is the case of models with two Higgs doublets (2HDM), the best known example 
being the minimal supersymmetric standard model in which two Higgs doublets are required to give mass to down-type 
quarks and charged leptons in one case, and up-type quarks in the other. In this framework, the ratios (jlj depend on 
the mass of the charged Higgs and the ratio /3 of the two Higgs doublet VEVs, and no choice of such parameters 
allows to simultaneously reproduce the experimental data on TZ{D) and TZ{D*) [T]. Variants of the 2HDM 01 [7], 
together with other models providing explicit flavour violation [3], might explain the measurements ([I]); however, an 
enhancement of the purely leptonic B decay rate is generally implied. 

In this paper we reconsider both the above mentioned issues. We reanalyse the SM prediction for B — > D''*^lvi^ 
specifying the main sources of uncertainties and possible improvements. Our results confirm that the most significant 
deviation is for TZ{D* ). Then, we scrutinize the effects of possible NP contributions in the effective weak Hamiltonian 
having a structure able to affect the ratios ([ij but leaving the pure B leptonic modes unchanged. In particular, we 
focus on a NP operator constructed from tensor quark and lepton currents. Such a kind of operators have been also 
investigated in [6j and [9j, but we devote the main attention to differential distributions, namely the lepton forward- 
backward differential asymmetries, in which the sensitivity to the new Dirac structure is maximal, as emphasized in 
[F for different operators. Although there are scenarios in which tensor operators are generated, in our analysis we 
do not rely on explicit models: our purpose is to identify physical observables having a mild sensitivity to hadronic 
uncertainties, which therefore can be used to unveil effects easier to interpret. It is only worth mentioning that 
these operators emerge, for example, in models with new coloured bosons carrying both lepton and baryon quantum 
number (referred to as leptoquarks, LQ): SU(5)gc/t I10]j Pati-Salam SU(4) |21j, composite [32], superstrings [23] and 
technicolor models [21]. In the most general formulation of these models scalar operators may also occur. Leptoquarks 
couple to quarks and leptons and, from limits on flavour changing neutral currents, preferably to those within the 
same SM generation. Searches for leptoquarks decaying to 2t and 26 jets, performed by the CMS Collaboration 
at the CERN LHC, bound (preliminarly) the mass of a possible scalar leptoquark to M{LQ) > 525 GeV, and to 
M{LQ) > 760 GeV for a vector leptoquark [2S]; other bounds can be found in [26] . 

In our analysis of semileptonic B decays, we first consider D and D* mesons in the final state, and then turn to 
the interesting case of final states with excited positive parity charmed mesons. 



II. EXCLUSIVE b clvi DECAYS 



We consider the h — ?> civi effective hamiltonian comprising the SM term and an additional operator [SI [S] : 



(4) 



Gf is the Fermi constant and Vcb the CKM matrix element, is the relative complex coupling of the new tensor 
term with respect to the SM one. It is assumed that the main coupling is to the heaviest lepton, hence we set = 
for £ = e, n and et = e-r- This coupling can be bound experimentally, so that the effects of the new operator can 
be scrutinized in physical observables which, in general, are expressed as a SM, a new physics and an interference 
contribution. For example, the differential B{p) — > Mc{p')£{pi)De{p2) decay rate, with Mc a charmed meson, reads: 



dT 



{B ^ M,£i?,) = C{q^ 



aq^ SM aq'^ 



dT 

NP dq^ 



{B MJi>i 



INT 



with q = p — p' and C{q^) defined as 



1 - 



(5) 



(6) 



A(a;, y, z) = + -\- — 2{xy + xz + yz) is the triangular function. To compute the three terms in ^ we need the 
relevant hadronic matrix elements. 



A. B -^ D£vi 



The hadronic matrix elements in B — >■ D£vi can be parametrized in a standard way. 



< D{p')\c-f^h\B{p) > = Fi{q^){p + p') 



[Foiq^) - Fiiq^)] q^ , 



(7) 
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niB+mD ' ' ruB + rriD 

(with Ft = Gt from the relation a^^^^ = ^Cfj^^ap <^°'^), so that the three terms in (|5|) read 



dq 
dt 
dq 
df 
d^ 



2 (B ^ DiDe) 



2 {B ^ DlDi) 



m% 1- 



^2 ^ ' 3m2 



2g2 



[FT(<z') + GT(g' 



(B ^ Df!?^) = -3i?e[eT 



(8) 

(9) 
(10) 

(11) 



In the infinite heavy quark mass hmit, formahzed by the heavy quark efltective theory (HQET), the form factors in 
(7]|8| can aU be related to the Isgur-Wise function ^ The result is known [351 [30]: expressing Fi{q^) and Fo(9^) 
in terms of two other form factors hj^iw) and h^{w): 



rrir, — m. 



Fi{q') 



W)-F,{q^)\ ^ 



1 



1 



[{niB + mD)h+{w) - [lUB - mD)h^{w)] 



(12) 
(13) 



2^mBmD 

and defining the meson momenta in terms of four-velocities, p ~ ttlbv and p' — mjjv' , with w — v ■ v' and q^ — 
rn^ + rnjj — 2raBraB)W, at the leading order in the heavy quark and as expansion one has 

h+{w)^aw) , h^{w)^Q , (14) 

with S,{w) the Isgur-Wise function. Also the form factors in Q are related to S,(w) at the same order expansion: 

-^-^""^ - ' '^'-^ (15) 



Friq ) = Griq ) = , i{w) 



y/uiBrnD 

At the next-to- leading order, corrections must be taken into account, which at first are needed for the study of the 
decay in SM. We elaborate a determination of the functions h^, and ^ based on a combination of experimental and 
theoretical information. The experimental input comes from the BaBar analysis of i? — > Dfj,D^ |31j . the differential 
rate of which, neglecting the lepton mass, reads: 



dT 



dw^ ' 487r3 



^^'^^'''m|r3(l + rf{w^ - lf/^[FD[w)f 



with 



Fd{w) = 



1 - r 



rriD 



and r — . Using the parametrization 

rriB 



F, 



in terms of the variable 



(w) = i^D(l)|l - 8p?z + (51/9? - 10)z2 - (252p? - 84)z3| 
_ y/^ITfT- %/2 



(16) 



(17) 



(18) 



(19) 



from the fit of the product G^''^''''{w) = FD{w)\Vcb\ the BaBar CoUaboration provides the parameters G'^''^°''(l) = 
F£i(l)|Vcb| and p\. The outcome of the fit is slightly diflterent for B~ or B^ modes; we consider for definiteness the 
case [S] Q 



G 



BaBar , 



1) = (44.9 ± 3.2 ± 1.6) 10^ 



pI = 1.29 ±0.14 ±0.05 



(20) 



1 The average between the charged and neutral B decay modes is quoted as G^"^'"'(l) = (42.3 ± 1.9 ± 1.4) lO'^ , pf = 1.20 ± 0.09 ± 0.04. 
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FIG. 2. Isgur-Wise function S,{w) (times |Vc6| x 10^) obtained using the BaBar data on B° -> D+rui (left) and the Belle data 
on — >■ D*^£~ue (right) . The width of the curves is due to the errors in the parameters fitted in the two cases and to the 
uncertainty on A and Qs in the determination of the form factor. 



This result can be translated into a determination of Ci^): expressing the form factors h±{w) in terms of the Isgur- 
Wise function and including the as and l/mb,c corrections worked out by M. Neubert in 29J and by I. Caprini et al., 
in 132,: 



h_{w) 



w + 1 



{C2 + C3) + (£6 + e,)Li 



w + 1 



(C2 - C3) + (ec - eb)^4 



(21) 
(22) 



with £6 = — , 

ZTOf, 



2toc 



The coefficients Ci.2,3 and Li are collected in appendix A Ci account for the perturbative 







corrections, Li for the heavy quark mass corrections and depend on the hadronic parameter A, the difference between 
the heavy meson (B, D) and the heavy quark (6, c) mass in the heavy quark limit. We use = 4.8 GeV and 
rric = 1-4 GeV and a conservative value A = 0.5 ± 0.2 GeV [29 , so that the uncertainty in A encompasses the error 
on A/to;, and A/ttIc. The Isgur-Wise function ^(w) resulting from 



G 



BaBa 



h+{w) - \^ h^{w) 



(23) 



is depicted in fig [2] (left panel). 

The form factors needed for analysis of the mode with r can be separately derived using again Eqs.( [2Tp2| ): 



\Vcb\h^{w) 



A{w) 



A{w) 



G 



BaBar 
BaBar / 



'iw) 



(24) 
(25) 



l+r 



with A = h-/h^. For the matrix elements of the tensor operator, we use ^(w) also in (15 1. In the standard model, 
the results for the semileptonic i?° — >■ branching fractions can be quoted as 



SM 



SM 



(2.15 ±0.45) X 10^2 
(0.70 ±0.12) X 10-2 



and, taking the correlation between the predictions for i and t into account, 

6(5° -> D+T'Vr 



SM B{B" D+£-De) 



0.324 ±0.022 



SM 



(26) 
(27) 

(28) 



The SM prediction for 7iP{D) deviates from the measurement ([T]) (with statistic and systematic uncertainties combined 
in quadrature) by about 1.5(t . The deviation is smaller in the charged TiT {D) case. 
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The stability of (28) against changes of the input information on form factors is noticeable: sensitivity to l/mg 
corrections can be estimated varying A, and this modifies the central value at a few per mille level. Sensitivity to the 
radiative corrections can be assessed changing the scale in as indicated in appendix |^ and also these corrections 
are not effective. Since the value at zero recoil G^'^^°''^{1) cancels out in the ratio, the main uncertainty in (28) comes 
from the error on the parameter experimentally determined. The value of TiP{D) coincides with the one obtained 
using the form factors Fi and Fq from lattice QCD with finite quark masses |6] . 



B. B ^ D'lvi 

While the results for TiP{D) and TZ^{D) do not display a statistically significant deviation from the SM expectation, 
the case of 72.° (D*), TZ~{D*) is quite different. The standard parameterization of the B ^ D* matrix element in 
terms of form factors is 



<D*{p',e)\-c^^{l-^5)b\B{p) > 



niB + mo 



"1r - ml,. 



(29) 



(with the condition Ao{0) = ^^^^i^^Ai(O) - ""^ "'^' ^2(0)) and 



2mD* 



< D*{p',e)\ccT^,{l - To)b\B{p) > = ro(q')- 



£ ■ q 

(mB + ruD* 



uaiiP"'p^ + ri(q^)e^^.Q/3p"e*'^ + T2{q'^)e^^ai3P°'e- 



niq'){e;p. - elp^) + T^^q'W^p', ~ elp'^) 
e*-q 



{Pt^P'u - PuP't,) 



with e the D* polarization vector. We choose the helicity basis for D 



(30) 



TO_D. V2 



(31) 



with E' andp' the D* energy and three- momentum in the B rest frame {E' = y^mjj, + jp'P and \p'\ — X{m'^Q , m^jj, , q'^)/2mB) 
The conditions ■ p' = and • e^,;, = —5ab, with a,b = L,±, are fulfilled. The differential decay rates for the 
longitudinal and the transverse D* polarization in terms of form factors are obtained from 



dTL 
dq 



2 {B ^ D*eDe) 



dTL 

dq 



2 {B ^ D*iDi) 



NP 



dTL 
dq 



2 (B ^ D*£Pe) 



INT 





6A( 






2g2 







,2 2 „2\2 



' q' 



I I \( 2 2 2\ \ I 2\ ^H'^B' "'D* J y ) A f 2\ 

(mB + mD')[mB - mB,. -q )Ai[q ) ■ ^2(17 ; 



A(m|,m|).,g2) 



2mj 
3me 



X{m%,ml,.,q^) 
mo'imB + mo')' 



rToiq^ 



rriB + mo* 



m%+ml,, -q^~ ,0 



(32) 



mo* 



Ti(g^) + 4mo.T2(g^) 



(33) 



4(m_B + mo- 



{mB + mo.) (mg ~ mjj, - q )Ai{q ) - \{mg,mjj,,q )A2iq ) 



X{m%,mjj,,q') ^,2n , 2(m|j + mf). ~g^) ~ 



mjj,{mB + mo') 



T,{q') + miq') 



(34) 



X{m%,m%,,q^) 
{mB + mD'Y 



[V{q')? 



(35) 
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NP 



(36) 



dg2 



INT 



I , 2 2 2n < / 2\ '^("^Bi '^i)* J 9^) Tr/ 2\ 

ms + mD')(mB - nijj, -q )Ai{q ^V{q ) 



[T,iq^)+Uq')] , 



to be multiplied by the factor C{q^) in ([6|. We have used the combinations 

Uq^)^To{q')-n(q^) 
fiiq')^T,{q')+Uq') 
f2iq')=T2{q')+niq') . 



(37) 



(38) 



At the leading order in the heavy quark expansion, the form factors in (p9|) and (30) are related to the Isgur-Wise 



function, while other contributions appear at the next-to-leading order. Analogously to the decay to D, one expresses 
V and Ai in terms of form factors hy and , 



V[1 ) = ^ , hv(w) 



2 77113. 
^i(g^) = ^mBniD* 



w + l 



iTiB + rriB' 



-hA^w) 



A2{q^) = 
^(9') = 



2^mBmD' 
1 

2y/mBmD' 



hA^iw) + —^hA^iw) 
niB 



msiw + l)hAi{w) — (niB — mo'w)hA2{w) — {vtlbW — mD*)hA3{w) 



(39) 



with q^ =^ w?n + m?). — 2mB'mD'W- Including as and — and — corrections, the relations have been worked out 

rrib rUc 



hviw) = [Ci -I- ec{L2 - Ls) + eb{Li - U)] Ciw) 
hAiM = 



' w + l 



w + l 



aw) 



hA.iw) ^ [Cl + e,{L, + Le)] ^w) 

hA,{w) = [Cf + Cl + e,(L2 -L3-L5 + Lq) + efc(Li - L4)] i{w) 



(40) 
(41) 

(42) 
(43) 



The expressions of C^, which incorporate the radiative corrections, and Li are collected in appendix [Xj the Li terms 
account for the ©(l/mg) corrections in the heavy quark expansion, and are determined from QCD sum rule analyses 
of the subleading form factors [29 . On the other hand, the relations of the form factors Ti in ( [30| to (,{w) in the 
heavy quark limit are: 



To(g')=T5(<z2) = 
Ti(g^)=T3(g^)= 

V 



aw) 



(44) 



TTlB 

ruD' 



aw) 



we use these expressions in the analysis of the tensor operator. 

Let us focus on the SM. Due to the heavy quark spin symmetry a unique form factor describes both _B — > I? 
and B — > D* transitions, so that we could use the Isgur-Wise function found in the previous section. To partially 
take into account the different experimental systematics, we choose to use the determination of ^ obtained by Belle 
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Collaboration from the analysis of B' 
mass, is 







D 



*+ 



dT 

dw 



(J5 ^ D*lvt) = 



GlW, 



], for which the differential decay rate, neglecting the lepton 



cb\ 



487r3 



{mB -mD> frn]j,g{w)T'^{w) , 



(45) 



with 



g{w)T^{w) 



(w) \Jw'^ - 1 {w + 1) 
(1 -r*)2 



, w — 1 
w + 1 



l + (l-i?2H) 



W — 1 



(46) 



In (461 r* = 



TUB 



, and C/, i?i and i?2 are given by 



Q{'w) — 



l(w + l)' 



1 — 2wr* 
(1 - r 



i?i(«;) = (i?*) 
i?2(zii) = (i?*) 



2^ + 1 V^(w) 



Ai(w) 



(47) 



2''« + 1 ^2(w) 



Ai(ii;) 



with R* = 2 ^ 
parametrization [52] 



The three unwnown fimctions in (46 



i?i(w) = i?4(l) - 0.12 (w - 1) + 0.05 (w - 1)2 
R2(w) = i?i(l) + 0.11 (w - 1) - 0.06 (w - 1)2 

(with z defined in (19l). The fit of the parameters in ( [48p0 1 is quoted as |3i 

J-(i)|v;^| ^ (34.6 ± 0.2 ± 1.0) x 10"^ 
p2 ^ 1.214 ±0.034 ±0.009 
i?i(l) = 1.401 ±0.034 ±0.018 
i?2(l) = 0.864 ± 0.024 ± 0.008 . 
From these expressions one can reconstruct 'f(w), 



47 1 have been determined by Belle adopting the 
(231p2 



91)z^ 



(48) 
(49) 
(50) 



(51) 



(52) 



with defined through Eq.(41 ). The fit provides us with the determination depicted in figl2^ (right panel). Through 
Eqs.( 40|42|43 ) the form factors hy, Ka^ and Ha^ can be reconstructed including the NLO ifniQ and Us corrections, 
and also B\B^ — )■ D*^t~vi) can be computed. The results are: 



BiB"" - 



> D*+t 
D*+T- 



SM 



SM 



(4.62 ±0.33) X 10-2 
(1.16 ±0.08) X 10^2 



and, taking the correlation between the predictions for the £ and r mode into account, 

_ S(50 ^ D*+T-Vr) 



SM 



B{B" D*+l'vi) 



= 0.250 ±0.003 



SM 



(53) 



(54) 



The result ( 54 ) deviates from the measurement in ([I]) (with statistic and systematic errors combined in quadrature) 
by 2.3cr. It coincides with the one in [2l[7l|9j, due to the stability of the ratio TlP{D*) against changes of the input 
parameters: varying the central value of A and of the quark masses by 30% produces less than 1% variation in the 
result. The radiative corrections, changing the scale in as mentioned in appendix|Xj do not produce an appreciable 
variation of the result. On the other hand, in the individual branching fractions there is a mild sensitivity to A: 
setting this parameter to zero (i.e. ignoring l/mg corrections) the branching fractions in (53 1 are reduced by about 
5%. In the charged case, there is a deviation of 1.8cr between the SM prediction for TZ{D*y&nA the measurement in 
0. 
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III. EFFECTS OF THE TENSOR OPERATOR ON H{D^*^) AND OTHER OBSERVABLES 

If the tensions in TliD) and TZ{D*) are due to NP effects, it is interesting to investigate the new operator in the 
effective Hamiltonian (Wj) which affects the observables in B — >■ D^*\vr transitions, focusing on the signatures with 
minimal dependence on hadronic quantities. As done in [2HZ1IS]) T^{D) and TZ{D*) data allow to constrain the values 
of the new effective dimensionless coupling. In our case ct is bounded as shown in figjS] Using the parameterization 

tT ^ larle'" + ero , (55) 

the tightest bound to e-ro and \aT\ is obtained from the measurement of TZ{D*), while the combination of Tl{D) 
and TZ{D*) data fixes the range of the phase 9. We select the overlap of the two regions determined by TZ{D) and 

TZ{D*) both at Icr. In this overlap region, the function x^i^r) = ( j^'^^AnJoy^p — ) + ( j^^^An{D')"'p^ — ) values 
running between 1.51 and 1.75. This permitted range of et is represented as 

i?e[eTo] =0.17 , /m[eTo] = , 

larl e [0.24, 0.27] (56) 
9 G [2.6, 3.7] rad 

and is also depicted in fig|3] Varying the effective coupling in this region we can analyze the impact of the new 
operator on various differential distributions. 

We start with the longitudinal and transverse D* polarization distributions in _B ^ D*TVr- We consider the decay to 

a D* with definite helicity, with differential decay width — -4— for the three cases L, ±. We define = H — , 

dq^ dq^ dq^ dq^ 

and show in fig|4]the differential branching fractions. The uncertainty in the distributions reflects the uncertainty on 
the parameters of the Belle Isgur-Wise function, on A and, in the case of NP, on e^. While the shape of the distributions 
are slightly modified from SM to NP, the maxima increase, a consequence of the increase of the branching fractions. 

The differential decay width distributions for D and D* (summed over the D* polarizations) have been measured 
by BaBar |34j, and can be compared to the SM and the NP scenario predictions. Once normalized to the total number 
of events, not only the SM distributions are compatible with data, as remarked in [SI] , but also the distributions in 
the considered NP scenario agree with measurements, as one can argue considering figj5j This confirms that the shape 
of such distributions does not allow at present to select between these possibilities, and other observables should be 
analyzed for a more efficient discrimination. 

Other observables are the longitudinal and transverse D* polarization distributions in _B — > D*TVr normalized to 
B ^ D*£Di. They are defined as 

rD* (2. _ dVL.TjB ^ D*T9r)ldq^ 




Re[er] Re[er] 



FIG. 3. (left) Regions in the {Re{eT), I'm{eT)) plane determined from the experimental data (to 1 and 2a) on TZ{D) (large 
rings) and TZ{D*) (small rings), (right) Region corresponding to values of between the minimum (indicated by the star), 
1.55 (yellow, light) and 1.65 (orange, gray) and 1.75 (brown, dark). 
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The SM predictions are shown in fig. [g] together with the modifications induced by the tensor operator. At large 
the observables are enhanced by 30 — 50 %, a noticeable effect. Furthermore, at odds with scenarios in which only 
Rl is affected by new physics [7], in the case of the tensor operator both the longitudinal and the transverse Rl and 
Rt distributions are modified. 

The longitudinal and transverse polarization fractions of the D* meson 



dVL,T{B ^ D*TVr) ( dT{B ^ D*TUr) 



dq^ 



-1 



(58) 



are shown in figjTj Both the SM and NP predictions are affected by a small error, since in the heavy quark limit 
the observables in ([58]) are free of hadronic uncertainties, due to the cancellation of the form factor S,{w) in the ratio. 
The residual uncertainty reflects that on A which controls the l/mg corrections. The uncertainty on A also enters 
in the curves obtained in the NP scenario in combination with ex- In SM, FL{q^) ranges between 0.75 at low q^ and 
about 0.35 at high squared momentum transfer; in NP in the allowed region of et, FL{q^) is between 0.35 and about 
0.65 at low (f' , while this observable converges to the SM value at high q^ . The SM predicts a dominant longitudinal 
polarization at small q^ , in NP the longitudinal and transverse polarizations have similar fractions up to q^ = GeV^ . 
An important observable is the forward-backward ApBiq^) asymmetry in B ^ DtVt and B — > D*tDt, defined as 



dT 



dV 
dq^ 



^ dq'^d cos 6 



(59) 



> 
O 

O 
X 

ta I 0.04 




> 

o 

O 
X 

13 1 0.05 




[GeV^] 



[GeV^] 



FIG. 4. Differential branching ratios with polarized D*: 



dB(B D*TVr)L , dB{B ^ D'rUr)T , . , ^> „, , 

^ — (left) and ^ — (ngnt). The lower 

dq^ dq^ 

(blue) bands are the SM prediction, the upper (orange) bands include NP effects. In SM, the uncertainties on the parameters 

of the Isgur-Wise function in Eq.(51|, together with the errors on A and as are included. In the case of the NP curves, the 



uncertainty on is also considered. 




q' [GeV'] 



q' [GeV'] 



FIG. 5. 



(left) and 



dq^ dq^ 
shaded histograms) compared to BaBar data (points) 



(right) distributions in the NP scenario (for the central value of et, 
HI; the distributions are normalized to the total number of events. 
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FIG. 6. D* polarization ratios 7?^ (g^) (left) and (q^) (right) defined in (571. Notations are the same as in fig|4] 



where 9i is the angle between the direction of the charged lepton and the Z)*-*-* meson in the lepton pair rest frame. 
We use the notation 



2, _ 1 3C(g2) J 



(60) 



with C{q^) defined in (l6| and the three terms in the parentheses given for D and D*: 



• D 



• D* 



Afiiq') = 8Fo(/)F,{q'){ml - m|,)^ (^1 - ^) X'/\ml,m^q') 



A^iq') - 

-^^^5^(9") = -&ReieT)Fo{q^)[FTiq') + Griq'MmB - mo 



Af^iq'^ 



)me^l-'^y'^'iml,m^,q^) (61) 



(1-^) X'/\ml,ml.,q') 



\mlAQ{q^) [A^{q^){mB + mD,f{ml - ml,, - q^) - A(m|, m^,. , g^)^^!'?^)] 
-AmD-{mB+mD^)q^A^{q^)V{q^)^ (62) 
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0.6 
0.5 
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FIG. 7. Polarization fractions Fl(q^) (left) and FT(g^) (right) for B D* tz/t defined in (|58k. Notations are the same as in 

figS 
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0.357 ±0.002 


-0.040 ± 0.003 


0.315 


0.026 


0.24 


0.07 


Afb 


0.40 ±0.005 


0.048 ±0.013 


0.30 ±0.005 


0.08 ±0.01 


0.21 ±0.003 


0.14 ±0.01 



TABLE I. Integrated forward-backward asymmetry for the considered decay modes. The first hne reports the SM resuhs, in 
the second hne the effect of the tensor operator is included. 



{ml - ml.){T,{q') + Uq')) + q^{T,{q') - Uq')) 



(63) 



Afb (r) = -4i?e(eT)™£ ( 1 " ^ 



\l/2/'™2 2 



mB + mo* 



A(m|,m|,,,g^) 
■niD* [mB + mD*)' 



rriD' 



niB + mD')Aiiq^){fiiq^) + faC?')) 
fi(92)+4mz3.f2((?2) 



g2(fi(g2) _ f2(g2)) ± (m| - ml,)if,{q^) ± f2(g2))] } 



(64) 



In fig|8]we plot ApBiq^) for B DrDr and i? D*tVt. The SM prediction is affected by almost no theoretical 
uncertainty, because of a nearly complete cancellation of the hadronic parameters in the ratio. In the case of NP, 
we have taken into account also the uncertainty on and \aT\- The SM curve lies in both cases below the NP 
distribution for all values of q^ . The most interesting deviation concerns the D* mode: the SM predicts a zero for 
Afb at q^ ~ 6.15 GeV^ in the NP case the zero is shifted towards larger values q^ G [8.1, 9.3] GeV^. Even though the 
experimental determination of the zero of the forward-backward asymmetry is challenging, this observable effectively 
discriminates SM from the NP model. The integrated asymmetries, obtained integrating separately the numerator 
and the denominator in ( 59 1 , are collected in Table |lj for D* , in the NP scenario the integrated asymmetry has 
opposite sign with respect to SM. 



IV. TENSOR OPERATOR IN B ^ D**tvi DECAYS 

The new operator in the effective hamiltonian Q affects other exclusive decay modes that are worth investigating. 
Of peculiar interest are the semileptonic B and Bg transitions into excited charmed mesons. The lightest multiplet of 
such hadrons, corresponding to the quark model p-wave — 1) mesons and gencrically denoted D*^*y comprises four 
positive parity states which, in the heavy quark limit, fill two doublets labeled by the (conserved) angular momentum 




[GeV^] [GeV^] 

FIG. 8. Forward-backward asymmetry Apsiq^) for B — )■ DrVr (left) and B — >■ D'tVt (right). The lower (blue) curves are 
the SM predictions, the upper (orange) bands the NP expectations. Uncertainty on A has been included and, in the case of 
NP, also on the parameters jarl and 9. 
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si = Sq + I (sq is spin of the light antiquark), hence S( = 1/2 or se = 3/2. The two mesons belonging to the first 
doublet, (-Dj^-jg, D'^s)i)^ ha.ve spin-parity — (0+,l+); the mesons in the second doublet have = (l+,2+) and 
are named £^(3)2)- members of the doublets, with and without strangeness, have been observed, and 

the two sf = 1/2+ states without strangeness are found to be broad, as expected |35| . 

In the heavy quark limit also the semileptonic B transitions to mesons belonging to the same charmed doublet can 
be described in terms of a single form factor. B decays to {Dq, D'l) are governed by a universal function denoted as 
'''i/2iw), B decays to (Di, D2) by the T3/2(w) form factor (the matrix elements are collected in appendix Ib|. There 
are several determinations of the Ti{w) parametrized in terms of the zero-recoil value Ti{l) (contrary to the Isgur-Wise 
function, Ti{w) are not normalized to unity a.t w = 1), of the slope and of the curvature q. In the ratios of branching 
fractions and asymmetries the zero-recoil value does not play any role, and this reduces the main dependence of the 
observables on the hadronic parameters. The present experimental situation needs to be settled, since the semileptonic 
B decay rates into {Dq, D[) exceed the predictions obtained using computed Ti{l); the origin of the discrepancy is 
still unknown, and could be related to the broad widths of the final charmed mesons, which determine a difhculty 
in the exclusive reconstruction, and to a possible pollution from other (e.g. radial) excited states. Semileptonic Bg 
decays to sf — 1/2+ cs mesons could clarify the issue, due to the narrow width of the strange charmed resonances 
[36) . On the other hand, the tensor operator produces precise correlations among various observables, therefore its 
effects could be distinguished from others. 

For definiteness, we use a QCD sum rule determination of t^/2{w) at leading order in as |37[ I38j. and of ti/2(w) at 



T-3/2(w) = T3/2(l) 1-pI/2{w-1) 



Tl/2{w) = ri/2(l) 1 - Pi/2(w - 1) + Cl/2(W - 1)' 



with 



2 

P3/2 



0.9 



r3/2(l) = 0.28 

ri/2(l) = 0.35 ±0.08 p^/2 = 2.5±1.0 Ci/2 = 3 ± 3 . 

The differential decay rates for B — > D**iDe can be written as in ([5|, see appendix [b| The ratios 



n(D* 



B{B 



BiB^D*£De) 



(65) 
(66) 

(67) 
(68) 



(69) 



and the analogous TZ{D[), TZ{Di) and TZ{D2) depend on the effective coupling er- This also happens in Bg — > Dl*£Dg 
transitions, in the SU{3) f symmetry limit for the form factors. 

In fig|9]for each meson doublet we show the correlation between the ratios ([69]) for B and Ss, together with the 
SM predictions {n{Dl) ,H{D'^)) = (0.077,0.100), {TI{DIq) ,n{D'^^)) = (0.107,0.112), {TZ{Di),TZ{D*) = (0.065,0.059) 
and {TZ{Dsi),TZ{D*2) = (0.060,0.055). The tensor operator produces a sizeable increase in the ratios TZ, which is 
correlated for the two members in each doublet. The hadronic uncertainty is mild: using the Ti functions in |40| . the 
results remain almost unchanged in the case of the = 3/2 doublet, while for si = 1/2 they are smaller by about 
25% in SM and in the NP case. The same effect is found using the form factors obtained by lattice QCD pTj . 

The differential forward-backward asymmetries in the case of the four positive parity charmed mesons are collected 
in fig 10 and the integrated ones in Table |l] While in B {D^, Di)TPr the forward-backward asymmetry does not 



discriminate between SN and NP, in the modes with D[ and D2 it is a sensitive observable: The inclusion of the 



tensor operator produces an enhancement of Afb with respect to SM for all values of q^. Moreover, in SM there is a 
zero which, in the case of _B D'^rDr moves towards larger values of q^, and disappears in B ^ D2TVr once NP is 
included. 

We close this section remarking that, while the tensor operator in Q does not affect the purely leptonic Be — T~h'r 
mode, it can have an impact on the transitions B^ — ^ {rjc^ J / ''P)t~ Vt and Af, — ^ h.cT~v^] therefore, sets of other 
observables can be identified and investigated, with precise correlated deviations from the SM predictions. 



CONCLUSIONS 



The detailed experimental information provided us on flavour physics shows an astonishing consistency with the 
SM predictions. The very few tensions identify possible paths to new physics searches. The BaBar anomalous 

B{B D^*'>Ti?r) 
B{B D(*)fiD^) 



enhancement of the ratios R{D'^*^) 



with respect to SM is one of these few cases. The analyses 
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R(0(.)i) 



FIG. 9. (left) Correlations between the ratios 71{D*^^^^) and 7?.(73j^)j) for mesons belonging to the {D*^^q, ^{a)i) doublet 
without (orange, dark) and with strangeness (green, light), (right) Correlation between 7?.(D(s)i) and Ti-{D*^-^2) for mesons in 
the (D(s)i, £^(3)2) doublet. The dots (triangles) correspond to the SM results for mesons without (with) strangeness. 



5^ 



0.5 
0.4 
0.3 
0.2 
0.1 
0.0 
-0.1 



> T Vr 



SM 



5 6 
[GeV^] 



[GeV^] 



5^ 



0.5 
0.4 
0.3 
0.2 
0.1 
0.0 
-0.1 



l^^D'X T Vr 




SM 



q' [GeV^] 




[GeV^] 



FIG. 10. Forward-backward asymmetry Afb for the decays B — > D^rVr (top, left), B — >■ D'^tVt (top, right), i? — >■ DiTUr 
(bottom, left) and B — ^ D2TVT bottom, (right) as function of . The solid (blue) curves are the SM predictions, the dotted 
(orange) bands the NP expectations. 



of R{D^*^) in specific models also evidentiatc the enhancement the pm^ely leptonic B — > tPt rate, for which data 
are better compatible with SM. A mechanism enhancing the semileptonic modes B — > D^*'>tDt with respect to 
B — >■ D^*^iJ,9ij_, leaving B — >■ tD^ unaffected, can be based on a tensor operator in the effective hamiltonian. We 
have bound the relative weight er of this operator and studied the impact on several observables, the most sensitive 
one being the forward-backward asymmetry in i? — >■ D*tVt with a shift in the position of its zero. If the anomaly 
in i? — J- D'^*^TPr is due to this NP effect, analogous deviations should be found in B to excited D transitions. The 
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ratios R for these mesons are enhanced with respect to SM, and the forward-backward asymmetry is a sensitive 
observable in the channels involving D[ and These signatures in exclusive semileptonic b — > ctDt modes make 
the understanding of the role of the new contribution to the effective weak hamiltonian feasible, a step towards 
possibly disclosing new interactions through flavour physics measurements. 
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Appendix A: Coefficients 



With the aim of providing the information useful to reconstruct the various B — ^ D^*") matrix elements, we collect 
here the expressions of the as and I/toq corrections in Eqs.( 21|22 1 and ( 40p3 ) worked out by M. Neubert and by I. 
Caprini et al. in [29l [32] . The functions Li (w) read as 

Li -0.72 (w- 1)A 
L2 ~ -0.16 (w - 1) A 
L3 ~ -0.24A 

2.4^0.24 A (Al) 

is --A 

3.24 



w + 1 



A 



The coefficients Ci are expressed in terms of Ci , 

CI 



4a< 



rf{w) 



a 



(5) 



a 



(5) 



2a, 
2a, 



Hi^) w, 



1 

'(5) ■ " 



with 



and 



(A2) 



rf{w) = 



log 



1 



(1 " logz„ =F 



2{w =F l)Zm{l ± 



l-2wzm + z"^ (1 - 2wz„ + 2;2J2 
- [(w ± 1) - 2w{2w ± l)z™ + [hw ± 211? T 1)^™ - 2z^]r/(w) . 
In ( A2|A4[ ) the lower signs refer to the index 5 (corresponding to the axial current). Ci reads: 



as{mc) 



a.hh{w) 



Zhh{w] 



ctsirric) 



log 



TTLr 



+ Zhh{w) + - [f{w) +rf{w)+ g{w)] 



with 



ahh{w) 
Zhh{w) 



wrf{w) — 1] 



27 

8 /94 



81 V 9 



(w- 1 

f{w) = wrfiw) - 2 - -4L=[L2(1 - wl) + {w'' ~ l)rj{w)] 
Vw^ - 1 

9{w) 



(A3) 



(A4) 



(A5) 

(A6) 
(A7) 
(A8) 



= [2.2(1 - ZmW_) - 2.2(1 - ZmW+)] - " I 2 'I [("^^ ~ + - 2:„) log(z„)] , 

1 \\- Zwz^ + z^) 



(A9) 
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and w± = w ± y/w'^ ~ 1. In the numerical analysis we set the scale /i = ^Jm^Pm^^ and investigate the sensitivity to 
higher order corrections varying this scale between /i/2 and 2/i. 

Appendix B: B — > D** matrix elements and differential semileptonic decay rates 

In the infinite heavy quark mass limit the B — > D** matrix elements can be defined in terms of two universal 
Txi2(w) and T3/2(w;) form factors: 

< i?*(p')|c7M(l - 75)^>|S(P) > = -2^mBmc.ri/2(^«) - v')^ ; (Bl) 
< ^o(p')|cCTp,.(l - 75)fo|S(p) > = l^niBmui ti/2(w) \_~t^,^afiV°-■v'^ ^ i(y^,v^ - VvV^^)\ ; (B2) 

< D\{^, e)|c7^(l - 75)^'|5b) > = -2^msmc;Ti/2M [-ie^^/S-e*"^^^'" " ("^ ^ 1)< + ' ; (B3) 

< ^'\{V ■, e)|ccr^i^(l - 75)^|^(p) > -l^mBmw^ T^jiiw) {-e^,.a^e*"(-i; - u')'' + i\(*^J<v - v')^ - e*(w - v)^,]} ; 

(B4) 

< e)|c7;.(l - 75)6|Sb) > - r3/2(^«) {*(1 + w)e^o,p,e*'^v^v'^ + {w^ - l)e; 

+ {e* ■ v) [3v^ - («; ~ 2^^] } ; (B5) 

< D^ip', e)\ca^,{l - To)b\B{p) > = r^/^^^) \ ~ {w ^ l)e^.a/3e*"(« + v'f + {e* ■ v)e^,^pv"v'^ 

v2 

+ + w) {etiv ~ - - - 3(e* • v){v^vl - v,v'^)] } ; (B6) 

< D;{p\ e)|c7^(l - -f5)b\B{p) > ^ ^m^m^ T3/2(w) { ~ ie^/j^^ (e*"'^fa) Vv'" 

+ i;^^ - (1 + i«) (e^^t,,) } ; (B7) 

< D*{p', e)|c(T^,(l - 75)^>|S(p) > = ^T^T^T^T^^V^Ty^iw) { - e^,^, (e*"''?;^) (^^ + 

+ * [5;:(^^ + - gliv + v'),] } ■ (B8) 

In the previous formulae we have set p — niB v, p' — raw v' and w — v ■ v' \ e is the polarization vector (tensor) of 
the spin 1 (spin 2) D** meson. 

The results for the SM, NP and interference contribution to the differential distributions in ([5| are given below 
for each of the four excited mesons. The relation between the squared momentum transfer and w \s, — + 
mjjrt — 2'mB'm£)**w, with rriB)-,* the mass of the charmed meson produced in the decay. The lepton mass has been 
taken into account, hence the formulae also hold for r. 

dV 



-^{B^D*eDe 



^^2v- • -u^^.. ^^^^AmBmD'jT,/2iw)]\w-l){q^ (l - ^) + + ^) [« + rnl^)^ - 2mBmD^] } 

— {B ^ D*£De) = 32|eT|'msmz54Ti/2(w;)]2(u.2 _ l) h + ^ (m| + m^,. - 2mBmD;w) (B9) 

dt 
dq 



„{B ^ D^tvi) = -A8Re{eT)TnBmD'[Ti/2iw)]'^{w'^ - l)mi{mB - mo*] 



INT 



• B^ D[£Di: 
dt 



dq 



.2 ' 



{B ^ D[£D,) = 4mBmD' [Tif2{w)Y{w - 1) 

SM ^ 
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NP 



dV 



dr 



INT 



dq 



2 {B ^ DiiDi] 



SM 



dT 
dq 



2 {B ^ DiiDe) 



NP 



dT 
dq 



2 (B ^ DiePi) 



INT 



• B^ D^epe 
dT 



dq 



2 (B -> D*eDe) 



SM 



dT 
dq 



2 {B ^ D;i9e) 



NP 



dT 
dq 



2 (B ^ 



/7VT 



{ 



(2w - 1) + 1 + 



2m^ 



(m^ + m^/ )3w — 2mBm£)j^ (2u;^ + 1) | 



32|eTpmsmi3. [ti/2(i«)]'(z« - 1) ( 1 + ^ j 

|(m^ + m^, )(5w — 1) — 2r7isTO£)/ [4 + — 

48i?e(eT)mBm£,/ [ti/2(w)]^(u) - l)mi[mBiw - 5) + m^,/ (5w - 1)] 



(BIO) 



mBniDi [T3/2i'w)f{w - 1)(1 + w)'^ 

. , / 2m? 
- 2) + 1 + 



(l - ^) - 2) + (^1 + [(m^ + m?,j3u; - 2msmDAy^' + 2)] } 

16|eTp™STOl3jT3/2H]2(u; - 1)(1 + wf h + ^ j 

{ [(m|+m^J(2u;-l)-2mBmDi(w'-w + l)] } (Bll) 
24:Re{eT)mBmDj^[T3/2{w)f {w - 1)(1 + w)^mi[mB{w - 2) + mD^lw - 1)] 



{g^ (^1 - ^ j (3w + 2) + (^1 + ^ j (m|+m|3.)5w-2mBmDj(3w2 + 2) } 
16|eT|'msmc* [r3/2(u;)]2(w; - 1)(1 + w)^ M + j 

I (m|+m^.)(l + 4w)-2mBmi5.(3 + w7 + u;2) | (B12) 
-24i?e(eT)mi3mi). [t3/2(w)]^(w - 1)(1 + w)2m£[ms(4 + w) - nio^ (1 + 4w)] 



The differential decay rates are obtained multiplying the above functions by the coefficient C{q^) in 
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